arXiv: 1501.04569v2 [hep-ph] 26 Apr 2016 


On the contribution 
of the double Drell—Yan process 
to WW and ZZ production at the LHC^ 


Mieczyslaw Witold Krasny® and Wieslaw Placzek^ 

Laboratoire de Physique Nucleaire et des Hautes Energies, 
Universite Pierre et Marie Curie Paris 6, Universite Paris Diderot Paris 1, 

CNRS-IN2P3, 

4 place Jussieu, 75252 Paris Cedex 05, France. 

^Marian Smoluchowski Institute of Physics, Jagiellonian University, 
ul. Lojasiewicza 11, SO-SfS Krakow, Poland. 


Abstract 

In this paper, we investigate consequences of an assumption that the discrepancy of 
the predicted and observed W'^W~ production cross sections at the LHC is caused 
by the missing contribution of the double Drell-Yan process (DDYP). Using our 
simple model of DDYP [Acta Phys. Pol. B 45, 71 (2014)], we show that inclusion 
of this production mechanism leads to a satisfactory, parameter-free description of 
the two-lepton mass distribution for 0-jet W~^W~ events and the four-lepton mass 
distribution for ZZ events. In such a scenario, the Higgs-boson contribution is no 
longer necessary to describe the data. An experimental programme to prove or 
falsify such an assumption is proposed. 
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1 Introduction 


The main motivation for our study presented in this paper was “the WW anomaly” at 
the LHC, he.deviations of the total cross sections for resonant W~^W~ production mea¬ 
sured by the ATLAS and CMS experiments from their theoretical predictions, 
as shown Table The cross sections measured by ATLAS at the collision energies of 
7 TeV and 8 TeV are higher respectively by factors of 1.16 and 1.22 than the theoretical 
predictions of the NLO QCD calculations obtained from the MCFM program [^. The 
corresponding factors for the CMS measurements are 1.11 and 1.22. 
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Table 1: Values of the total cross section for resonant W~^W~ production at the LHC measured by the 
ATLAS and CMS experiments and predicted by the theoretical NLO QCD calculations from the MCFM 
program. 


Although each of these deviations does not exceed 2cr, the fact that they are present 
in four measurements and all exhibit excess of data with respect to theory may indicate 
that there exist some extra processes contributing to the measured cross sections that 
have not been taken into account in the theoretical predictions. Such processes may have 
an important impact on the significance of the Higgs-boson signals at the LHC. 

Recently, the NNLO calculations for the W~^W~ production have been published [^. 
They predict increase of the theoretical cross sections at these energies by about 9% with 
respect to the NLO results, so they get closer to the experimental measurements but do 
not remove the differences completely. These calculations have been done, however, for the 
total cross sections only and not for distributions considered by the LHC collaborations in 
the Higgs searches. The question if the increase of the cross section is in the Higgs-signal 
region or in the Higgs-monitoring region, or both, which is critical to our analysis, remains 
thus open. Therefore, we shall not use them in this paper. We can do this in the future 
when the ATLAS and CMS collaborations apply them in their data analyses. 

Since the observed W~^W~ cross section discrepancies at 8 TeV are by a factor of 
~ 3 higher than the expected contribution coming from the Standard Model (SM) Higgs- 
bosons decays, special measurement procedures were used by the ATLAS and CMS col¬ 
laborations, see e.g. Refs. [^[^, to increase their sensitivities to the Higgs signal. The 
normalisation of the predicted background was rescaled to fit the data in the monitoring 
region, where the Higgs contribution is negligible. Then, the rescaled background was 
used in the Higgs-search region to determine the strength of the Higgs signal. 

In the case of the ATLAS experiment, the normalisation of the theoretical predictions 
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for the background was multiplied by a factor of 1.22 for the 8 TeV analysis of the 0-jet 
e/x events [^, where the sensitivity to the Higgs boson decays is the highest. This value 
is compatible with the discrepancy of the measured and predicted total cross sections for 
the resonant W~^IV~ production at 8 TeV, see Table [Tj The actual value of the rescaling 
factor for the CMS analysis is not explicitly quoted in the their papers. Therefore, the 
CMS data will not be used in the presented analysi^ 

In this paper, we shall make a bold assumption that the missing process accounting for 
the cross-section discrepancy, not considered so far in the calculations of the theoretical 
cross sections, is the double Drell-Yan process (DDYP) resulting from double-parton 
scattering (DPS) in proton-proton collisions. It should be stressed that the strength of 
DDYP can, at present, be neither directly constrained by experimental data nor predicted 
theoretically. 

If one includes DDYP as a contributor to the W~^W~ production processes, it is bound 
to contribute as well to the ZZ production processes with a fully constrained strength. 
The immediate question which may be asked is if after adding the DDYP contribution 
to the Higgs-boson searches background, both in the W~^W~ and ZZ channels, there 
would still be a need to include the contribution coming from the Higgs-boson decays 
or, putting it alternatively, to which extent the DDYP contribution, on top of curing 
the WW anomaly, could mimic the Higgs-boson signal in the 125 GeV Higgs-sensitive 
phase-space regions. 

To answer this question, we present a coherent analysis of the DDYP contribution to 
the H —)■ W~^W~ and H ^ ZZ* decay-channel backgrounds. We use the data corre¬ 
sponding to the total integrated luminosity collected so far at the LHC. Once the overall 
normalisation of the DDYP contribution is hxed to explain the WW anomaly, our DDYP 
model does not have any more free parameters, thus can be easily falsihed by comparing 
its predictions to experimental data. 

The DDYP contribution to the Higgs searches background in the H ^ ZZ* decay 
channel was already analysed in our previous work 10 . Using a simplihed model of DDYP, 


we have demonstrated the appearance of a peak in the four-lepton invariant mass, 
distribution in the ~ 125 GeV Higgs-signal region. This “Higgs-like” peak is not driven 
by the details of the DDYP model. It is generated by the interplay of a steeply falling mu 
distribution and the kinematical threshold effect driven by the experimental cuts on the 
outgoing leptons variables: the minimal transverse momenta of leptons and the minimal 
invariant mass of the opposite charge lepton pairs. These cuts are similar in the ATLAS 
and CMS analyses, and therefore result in similar DDYP peak positions within a 2 GeV 
interval. 

It has to be stressed that the claim of the ATLAS and CMS collaborations that 


more general comment is obligatory for discerning unavoidable caveats of the analysis presented in 
this paper. The LHC collaborations publish very rarely their detector-effect unfolded distributions. The 
Higgs papers are no exceptions here. This preempts a fully irrefutable justification of any external analysis 
of these distributions, including the one presented in this paper. All we can do is to try to minimise the 
impact of the necessary underlying assumptions on the event selection efficiencies and detector smearing 
effects. Their remaining impact can be evaluated only by the relevant collaborations. 
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DDYP can be neglected as a potentially alarming sonrce of backgronnd was based on 
the assnmption of nncorrelated: (1) longitndinal momentnm, (2) transverse position, (3) 
flavonr, (4) charge and (5) spin of the partons taking part in DDYP, and on the assnmption 
of the process-independent valne of aes, governing the strength of the DPS processes. The 
above, in onr view nnjustihed, assnmptions lead to a signihcant nnderestimation of the 
contribntion of DDYP to the Higgs searches backgronnd. As we argued in Ref. 10 , its 


contribution must be, given the lack of the adequate theoretical calculations, determined 
experimentally, e.g. by using the experimental methods proposed therein and in Ref. [^. 
There are three main reasons for writing this paper: 


1. In Ref. 10 , we applied our DDYP model only to the ATLAS data for the ZZ 


channel, while in this paper, we apply it to the ATLAS data for the WW channel, 
and after htting the normalisation factor, we use it, parameter-free, in the ZZ 
channel. 


The analysis in Ref. 10 was done only for the partial ATLAS data available at that 
time, here we use the full data collected by ATLAS in the LHC Run 1 to check if 
our DDYP model can still describe these data. 


3. We propose new experimental methods of testing DDYP at the LHC, see Section]^ 


The paper is organised as follows. In the next section, we present numerical results of 
our analysis. Section includes a discussion of the results as well as a general discussion 
of the interplay between the DDYP and gluon-gluon scattering contributions. It contains 
a proposal of measurements capable to elucidate the role of DDYP in WW and ZZ 
production processes. Finally, Section concludes our paper. 


2 Results 


The numerical results presented below have been obtained using the Monte Carlo event 
generator WINHAC 11 13 with the same model of DDYP and the same input parameters 
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as m our previous paper 

The starting point to our quantitative studies in this paper is a shape analysis of 
the two-lepton mass, mu, distribution and its £t to the -|- 0-jets” hnal-state data 

presented by the ATLAS Collaboration in a wide mjj range for the full data statistics at 
8 TeV in Ref. 

In Fig. (a), we present the results of the ATLAS Collaboration of Ref. for the 
collision energy of 8 TeV. The background and the Higgs decays contributions are shown 
separately. This plot reflects the necessity of including the Higgs contribution to obtain 
a satisfactory description of the mu distribution in the region of its small values, where 
the sensitivity to the 125 GeV Higgs boson is the highest. However, we would like to 
stress again that the original background prediction was, in this case, rescaled by the 
ATLAS Collaboration by the factor of 1.22 in order to get a good description of data in 
the monitoring region. 
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Figure 1: The mu distributions for the H — >• W^W~ searches at ^/s = 8 TeV: (a) the rescaled (ATLAS) 
background without (grey histogram) and with (solid/red line) the SM Higgs boson contribution compared 
with the ATLAS data (black dots) , and (b) the canonical background without (grey histogram) and 
with (solid/blue line) the DDYP contribution compared with the same ATLAS data (black dots). 


For the plot presented in Fig. [^(b), we re-normalise back the background distribution 
to its original, canonical theoretical predictions by dividing the background shown in 
Fig. a (a) by a factor of 1.22, and subsequently add the DDYP contribution with a 
normalisation of its prediction determined by minimisation of the x^/d-o.f. in the £t of 
the sum of the DDYP and theoretical background contributions to the data. The Higgs 
contribution is no longer necessary to obtain a satisfactory description of the data by the 
canonical background model if the DDYP contribution is included. 

Indeed, the values of x^/d-o.f. corresponding to the mu distributions in Fig. are 
1.2 for the “rescaled background plus the SM Higgs” scenario (with the p-value equal to 
0.26) and 0.8 for the “canonical background plus the DDYP” scenario (with p-value equal 
to 0.70). Both descriptions of the data are acceptable on a statistical basis, although 
the obtained values of the likelihood test give some better preference to the “DDYP plus 
canonical background” predictions. We stress again that in the latter case, the theoretical 
predictions of the SM background do not need to be rescaled to describe the data, as it 
was done in Refs. (^[^. 

Having determined the normalisation of DDYP cross section using the W~^W~ chan¬ 
nel, we can now study the corresponding DDYP contribution to the m^i distribution for 
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Figure 2: The m 4 i distributions for the H —>■ ZZ searches at = 8TeV: (a) without (grey his¬ 
togram) and with (solid/red line) the SM Higgs-boson contribution compared with the ATLAS data 
(black dots) [^, and (b) the (ATLAS) background plus the DDYP contribution (solid/blue line) com¬ 
pared with the same ATLAS data (black dots). 


the ZZ channel. In our model, the dominant contribution to the DDYP cross section 
comes form the qq excitations of the proton sea. Thus, the relative normalisation of the 
DDYP contributions for the W^W~ and ZZ channels is driven, on the modelling side, 
only by the relative strength of the uu and dd excitations, reflecting the ratio of the u 
and d PDFs (for more details, see Ref. |^). On the experimental side, it reflects the 
relative efficiencies of selection of W^W~ and ZZ events in their kinematical acceptance 
region^ The DDYP model predictions for the ZZ channel are thus fully constrained. 

Again, as for the W^W~ channel, the ATLAS data at 8 TeV 14 are presented hrst 
in Fig. (a) with the SM Higgs contribution and the canonical background contributions 
shown separately. In Fig. (b), we replace the Higgs contribution by the parameter-free 
DDYP model predictions. The values of y^/d.o.f. corresponding to mu distributions in 
Fig. I are 0.65 for the background plus SM Higgs scenario (with the p-value equal to 
0.94) and 0.95 for the background plus our DDYP predictions (with p-value equal to 
0.55). Both descriptions of data are statistically acceptable, however in this case, the SM 


^For the analysis of the W~^W~ and ZZ spectra, we select events according to the selection criteria 
used for the Higgs searches in the corresponding decay channels. We assume the same efficiency of the 
W~^W~ and ZZ event selection and neglect the detector smearing effects. 
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Higgs prediction is slightly more preferred. 

In contrast to the SM Higgs scenario, DDYP contributes not only to the m^i range 
near 125 GeV but also above the threshold of the on-shell ZZ production. We have 
calculated that with W^W~-f\x.ed normalisation, DDYP will result in 31 extra events in 
the range of 190 GeV < m^i < 250 GeV, which agrees at the la level with the ATLAS 
data excess of 20 ± 13.5 events over the SM background 14 . GMS also observes an excess 


of the data with respect to the background in this kinematical region at a comparable 
level of ~ 13% 
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3 Discussion 


Our numerical results presented in the previous section show that once the DDYP contri¬ 
bution is normalised to account for the missing contribution to the W^W~ cross section, 
one obtains a satisfactory predictions for the excess of events observed by the ATLAS ex¬ 
periment both in the mu and 7714 ^ distributions. These excesses of events were attributed 
by the ATLAS and GMS collaborations to the SM Higgs-boson decays to ZZ and W'^W~ 
with mn ^ 125 GeV. 

The presented results correspond to the two extreme scenarios explaining the excess of 
the experimental data over the SM background: the “SM-Higgs-only” and the “DDYP- 
only”. There is, of course, a possibility that the two could contribute together. For 
example, in the case of the distribution, the SM Higgs contribution could better 
describe the peak width in the data in the 125 GeV region, while DDYP can add some 
events in this region but also in the higher mass range. This would lead to a different 
values of the Higgs couplings to the W and Z bosons with respect to those quoted in the 
canonical ATLAS and GMS analyses which neglect the DDYP contribution altogether. 

The intriguing question is if the DDYP contribution could be sufficiently large to 
mimic fully the Higgs signal at the LHG. An equally important question is why it was not 
considered as important background source in the data analyses by the ATLAS and GMS 
collaborations. Let us start with the answer to the latter question before addressing the 
former one. 

The DDYP contributions were studied by both collaborations with the help of the 
standard parton shower Monte Garlo generators and were found to be negligible. The 
reason why these generators predict such small cross sections for DDYP is related to the 
assumption of a completely uncorrelated two single-parton scattering (SPS) modelling of 
two Drell-Yan processes in a single proton-proton collision. The canonical DPS models 
assume the value of the so-called effective double-parton scattering cross section: Ueff ~ 
15 mb. This parameter normalises the DPS cross section with respect to the product of 
two SPS cross sections, see e.g. 
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The above value has been measured, but only for the cases of DPS involving at least 
one gluon in a pair of colliding partons. There are at least three reasons why the value 
of (Teff may be signihcantly lower for the same flavour, opposite-charge and spin quark- 
antiquark pairs relevant to the ZZ and W^W~ production processes: 
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• The quark-antiquark excitations of the proton involves partons of the same flavour. 
For example, if an s-quark takes part in the production of one of the vector boson, 
its s partner is already present in the wave function of the colliding proton. This 
enhances by a large factor, with respect to the canonical picture in which the PDFs 
of both partons are considered as uncorrelated, the probability that it may take part 
in the production of the second vector boson. 

• The transverse-plane correlation length for the same flavour and opposite charge qq 
pairs is signihcantly smaller than for the gluon-gluon pairs because of the dominance 
of the local charge and flavour conserving proton excitations — note that at the 
hardness scale of the vector-meson pair production processes and for the x-region 
relevant to the Higgs searches, there is less than one hxed-flavour gg-pair excitation 
per incoming proton, while there is more than 100 gluon pairs covering uniformly, 
glueing together, the proton volume. No doubt, the strength of the gluonic DPS 
must reflect the transverse size of the proton, while DDYP should reflect rather the 
typical size of a quark-antiquark dipole. 

• The long-lived qq excitations within the proton which do not lead to the signihcant 
effects modifying its overall spin involve quarks and antiquarks of opposite spins. 
Collisions of such pairs produce always vector-boson pairs with the total spin equal 
to zero, as in the case of the Higgs-boson decays. 

All the above amplihcation effects cannot be calculated within the pres-ently available 
QCD perturbative methods. This, however, does not mean that they do not exist. The 
LHC provides the unique opportunity to measure the respective quark-antiquark flavour, 
transverse plane, longitudinal momentum and spin correlations in the proton. The XY- 
pair production processes, where X,Y G {■j*, Z,W, J/^,T}, provide an excellent exper¬ 
imental testing ground for the quark-antiquark correlation models, in particular if: (1) 
data are collected at two or more collider energy settings (allowing to separate experi¬ 
mentally the quark and gluon originating processes D). (2) the data are collected not 
only in the pp but also in the pA and AA collision modes (to control the transverse-plane 
parton correlations) and, most importantly, (3) the DDYP effects are measured both in 
the Higgs-signal and monitoring regions. In addition, the measurement of the relative 
strengths of these processes provides a clear experimental test of the robustness of the 
SM Higgs interpretation of the data with respect to alternative mechanisms of the elec- 
troweak symmetry breaking. As long as such experimental studies are not made, the 
DDYP model should be, in our view, considered as equally plausible as the Higgs model 
in explaining the source of the excesses of events in the ZZ and W~^W~ channels observed 
at the LHC. 

In this paper, we avoid giving any value of cxefr for our DDYP, since, in general, DPS 
does not simply factorise into the product of two SPS processes, especially when the DPS 
contribution is sizeable, as it is in our case; examples of that are shown e.g. in Ref. pT] . 
Even if it were possible, the ATLAS and CMS experiments do not provide the necessary 
information (detector efficiencies and smearing effects) to translate our normalisation 
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factor to (Jeff. At the generator level, this normalisation factor is of the same order of 
magnitude as in Ref. 
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Finally, let us address a more subtle question of the interplay between DDYP and 
the electroweak-boson pair production in the gluon-gluon fusion process. The latter pro¬ 
cesses are included in the MCFM program as well as in the gg2WW and gg2ZZ [I^ 
generators which are used by the ATLAS and CMS collaborations for theoretical pre¬ 
dictions of the respective SM background to the Higgs signal in its IVW and ZZ decay 
channels. These calculations include also the so-called crossed-box contribution in which 
two incoming gluons split into quark-antiquark pairs and then a quark (antiquark) of one 
pair interacts with an antiquark (quark) of the second pair, leading to DDYP. One may, 
therefore, think that the processes we consider in this paper are already included in the 
theoretical predictions of the SM background to the respective Higgs signals. 

In our opinion, the gluon-gluon fusion calculations include only partially the DDYP 
contribution and may even underestimate the included crossed-box diagram effects. Firstly, 
they take into account only on-shell incoming gluons which are purely left- or right-handed, 
while, as we argue in Ref. 10 , one should consider all possible gluon polarisations. For 


the spin-zero vector boson pairs mimicking the Higgs signal, a particular care must be 
given to the full set of processes producing on-shell and off-shell qq pairs with compensat¬ 
ing spins. Such pairs may be produced by longitudinally polarised gluons but also by the 
higher-twist effects. To our best knowledge these effects are not included in the existing 
theoretical calculations. 

As was shown in a detail in Ref. 19 , the crossed-box contributions exhibit a collinear 
singularity when the spins of the incoming on-shell gluons sum to zero. The typical 
collinear singularity is, within the leading-twist approach, damped in this case to the 
logarithmic (integrable) singularity due to the vector structure of QCD for the on-shell 
initial-state gluons, he.when they are purely left- or right-handed 19 . Off-shellness of 


longitudinally polarised gluons should, to some extent, damp the singularity, but will this 
reduce the enhancement due to the collinear quark-antiquark pair emission? 

The question which remains to be answered is not only how frequently the incoming 
gluons are longitudinally polarised but also what is the probability for the qq pair prop¬ 
agating over large distances, before annihilating into a vector-boson pair, to become a 
spin-zero pair due to soft colour interactions with the medium. In both cases, the total 
spin of the colliding qq dipoles is zero, leading to a spin zero configuration of the hnal 
vector-boson pairs. 

Another problematic issue is the question of the renormalisation and factorisation 
scales. In typical fixed order QCD calculations, as given in Refs. [5,17-19 , these scales are 
set equal to each other and taken as a hard-process scale. In the case of the processes under 
consideration, this scale is of the order of ~ 100 GeV. Is using such a high scale justified 
for processes in which incoming gluons split into almost collinear quark-antiquark pairs? 
In our opinion it is not. In such cases, a better choice for the argument of the running 
QCD coupling may be not the hard process scale but rather the transverse momentum 
of the emitted quarks, see e.g. Ref. [^. Generally, this kind of a scale is used in popular 
QCD parton shower generators. This is, however, often not the case in the hxed-order 










QCD calculations. The value of the running ag between the scales of ~ 100 GeV and 
~ 1 GeV is increased by a factor of ~ 5. Since the considered process involves Og, we 
can easily get the enhancement factor of ~ 25. The factorisation scale for such collinear 
splittings should also be set to a similar value. 

The effects discussed above should lead to enhancements in low px regions of pro¬ 
duced electroweak bosons. This is hard to observe in the case of the production, 

as the transverse momenta of neutrinos from the leptonic fV-boson decays cannot be mea¬ 
sured individually, but could be seen rather easily in the ZZ jZ^* processes through their 
charged lepton decay channels. Therefore, the latter processes can provide an important 
test of the interplay between the DDYP and gluon-gluon scattering contributions. 

4 Conclusions 

In this paper, we have investigated the consequences of the hypothesis that the double 
Drell-Yan process (DDYP) accounts for the excess of the measured W~^W~ cross section 
with respect to its theoretically predicted value. This assumption determines the absolute 
normalisation of the predictions of our simple DDYP model introduced in Ref. (^. This 
normalisation factor is the only free parameter of the model. 

We have demonstrated that adding the above absolutely normalised DDYP contri¬ 
bution to the canonical SM background is sufficient for a satisfactory description of the 
17141 spectra for the ZZ hnal state and 17121 spectra for the W~^W~ hnal state, with a 
comparable £t quality as in the model assuming the existence of the SM Higgs boson. 

We have argued that the DDYP contributions may indeed be signihcantly larger than 
that expected in the naive canonical DPS models because of the strong charge, flavour, 
longitudinal momentum, transverse position and spin correlations of the quark-antiquark 
pairs participating in DDYP. Such a possibility has not been excluded so far, neither by the 
theoretical calculations nor by the experimental measurements. We have presented some 
theoretical arguments in favour of such a DDYP contribution and proposed a measurement 
programme to test it experimentally at the LHG. 
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